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ABSTRACT 

Aims. We present FUSE spectra of the giant Hn region NGC604 in the spiral galaxy M33. Chemical abundances are tentatively derived from 
far-UV absorption lines and compared to those derived from optical emission lines. 

Methods. Absorption lines from neutral hydrogen and heavy elements were observed against the continuum provided by the young massive 
stars embedded in the H n region. We derived the column densities of H i, N i, O i, Si n, P n, Ar i, and Fe n, fitting the line profiles with either 
a single component or several components. We used CLOUDY to correct for contamination from the ionized gas. Archival HST/STIS spectra 
across NGC604 allowed us to investigate how inhomogeneities affect the final Hi column density. 

Results. Kinematics show that the neutral gas is physically related to the H n region. The STIS spectra reveal H i column density fluctuations 
up to 1 dex across NCG604. Nevertheless, we find that the Hi column density determined from the global STIS spectrum does not differ 
significantly from the average over the individual sightlines. Our net results using the column densities derived with FUSE, assuming a single 
component, show that N, O, Si, and Ar are apparently underabundant in the neutral phase by a factor of ~ 10 or more with respect to the ionized 
phase, while Fe is the same. However, we discuss the possibility that the absorption lines are made of individual unresolved components, and 
find that only P n, Ar i, and Fe n lines should not be affected by the presence of hidden saturated components, while N i, O i, and Si n might be 
much more affected. 

Conclusions. If N, O, and Si are actually underabundant in the neutral gas of NGC604 with respect to the ionized gas, this would confirm 
earlier results obtained for the blue compact dwarfs, and their interpretations. However, a deeper analysis focused on P, Ar, and Fe mitigates 
the above conclusion and indicates that the neutral gas and ionized gas could have similar abundances. 
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1. Introduction 

The interstellar medium (ISM) is mainly enriched by heavy 
elements produced by the young massive stars during many 
starburst episodes taking place over the star-formation history 
of the galaxy. The immediate fate of metals released by these 
massive stars in the H n regions where stars recently formed, 
has not yet been settled. Kunth & Sargent (1986) have sug- 
gested that the H n regions of the blue compact dwarf galaxy 
IZw 18 enrich themselves with metals expelled by supernovae 
and stellar winds during the timescale of a starburst episode 
(i.e., a few 10 6 yr). However, observational evidence (see e.g., 
Martin et al. 2002) shows that metals might be contained in a 
hot phase reaching the halo of galaxies, before they could cool 
down and eventually mix into the ISM. The issue of the pos- 
sible self-enrichment of H n regions is essential to understand- 
ing the chemical evolution of galaxies, since H n region abun- 
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dances derived from the optical emission lines of the ionized 
gas are extensively used to estimate the metallicity of galax- 
ies. If there is self-pollution of these regions, the derived abun- 
dances would no longer reflect the actual abundances of the 
ISM. 



One approach to studying possible self-enrichment and, 
more generally, the mixing of heavy elements in the ISM is to 
compare the abundances of the ionized gas to those of the sur- 
rounding neutral gas. A first attempt was made by Kunth et al. 
(1994) using the GHRS onboard the Hubble Space Telescope 
(HST). These authors derived the neutral oxygen abundance in 
IZw 18 using the O i A 1302 line arising in absorption and using 
the H i column density and velocity dispersion from radio ob- 
servations. They found an oxygen abundance that is 30 times 
lower than in the ionized gas. However, the result remains in- 
conclusive since the O i line is likely to be heavily saturated, so 
that its profile could be reproduced with a solar metallicity by 
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choosing a different dispersion velocity parameter, as pointed 
out by Pettini & Lipman (1994). 

The Far Ultraviolet Spectroscopic Explorer (FUSE; Moos 
et al. 2000) gives access to many transitions of species arising 
in the neutral gas, including Hi, Ni, Oi, Sin, Pn, An, and 
Fe n, with a wide range of oscillator strengths for most of the 
species. Hence it becomes possible to determine abundances in 
the neutral gas with better accuracy. 

A surprise came from recent FUSE studies of five gas-rich 
blue compact metal-poor galaxies (BCDs): IZw 18 (Aloisi etal. 
2003; Lecavelier et al. 2004), Markarian 59 (Thuan et al. 2002), 
IZw 36 (Lebouteiller et al. 2004), SBS 0335-052 (Thuan et al. 
2005), and NGC625 (Cannon et al. 2004). In these galaxies, ni- 
trogen seems to be systematically underabundant in the neutral 
gas with respect to the ionized gas of the H n regions. The oxy- 
gen abundance seems either lower in the neutral gas or similar 
to the one in the ionized gas. This picture, if true, could offer 
a new view of the chemical evolution of the ISM in a galaxy, 
especially for the metal dispersion and mixing timescales. It 
could, however, suffer from uncertainties such as ionization 
corrections, depletion effects, or systematic errors due to the 
multiple sightlines toward individual stars and, possibly, multi- 
ple H ii regions within the spectrograph entrance aperture. 

In this respect, nearby giant H n regions in spiral galaxies 
provide an interesting case, since only one H n region fits the 
aperture. Moreover, their resolved young stellar population can 
be investigated to analyze the stellar continuum and in order to 
understand and correct for the effects of having multiple sight- 
lines contributing to the absorption line profiles. NGC604 is 
the first nearby giant H n region we have investigated for this 
purpose. It is a young star-forming region in the dwarf spiral 
galaxy M33 with an age of 3-5 Myr (e.g., D'Odorico & Rosa 
198 1 ; Wilson & Matthews 1995; Pellerin 2006). NGC604 is the 
brightest extragalactic Hn region in the sky after 30 Doradus. 
It is located about 12' from the center of M33. At a distance of 
840 kpc (Freedman et al. 1991), 1" corresponds to 4.1 pc. The 
nebula has a core-halo structure, where the core has an opti- 
cal diameter of ~ 220 pc and the halo a diameter of ~ 440 pc 
(Melnick 1980). The optical extinction in NGC604 appears 
to be correlated with the brightness. It is highest toward the 
brightest regions (Ay = 1.7-2.8 mag), while the average ex- 
tinction over the nebula is Ay ~ 0.5 mag, suggesting that dust 
is correlated with the ionized gas (Churchwell & Goss 1999, 
Viallefond et al. 1992). 

We first describe the observations in Sect. [2] and the data 
analysis in Sect. [3] with particular emphasis on the influence 
of the source extent on the absorption line profiles. In Sect. |4] 
we analyze the diffuse molecular hydrogen content, while in 
Sect. |5] we infer the neutral hydrogen column density. Metal 
column densities in the neutral gas are discussed in Sect. [6] 
We estimate the neutral gas chemical composition in Sects. 
and|8]and eventually compare it with the ionized gas of the H n 
region. Finally, in Sect. |9] we investigate how the presence of 
hidden saturated unresolved absorption components can be a 
major problem for determining the column density. 




Fig. 1. FUSE and HST/STIS apertures are plotted over the optical im- 
age from the POSS2 survey (B band). We superimposed the most mas- 
sive stars in black as observed with STIS at 2000 A. The arrows show 
the dispersion direction for each slit. North is up. 

Table 1. Log of the observations. 





FUSE 


FUSE 


HST/S77S 


IUE 


PID 


B018 


A086 


9096 




Date 


12/2003 


09/2001 


08/1998 


1979-1984 


Range (A) 


900-1200 


900-1200 


1150-1730 


1150-1950 


Exp. (ksec) 


7 


13 


2 


4.8-22.8 


Aperture 


4"x20" 


30"x30" 


52"x2" 


«10"x20" 




(MDRS) 


(LWRS) 


G140L 


SWP 


Resolution" 


0.07 A 


0.07 A 


2 A 


5 A 



a Spectral resolution for a point-like source. 



2. Observations 

The ionizing cluster of NGC604 was observed with the FUSE 
telescope through the MDRS and LWRS entrances (see the log 
of the observations in Table 0. The source as observed in the 
far-UV, with an apparent size of ~10"xl5" (see Fig.0, should 
be considered as extended relative to the size of the aper- 
tures. Data were recorded through the LiF and SiC channels 
(~ 1000-1200 A and ~ 900-1 100 A resp.). These channels are 
independently calibrated and provide redundant data, making 
it possible to identify possible instrumental artifacts. The data 
were processed with the CalFuse 2.4 pipeline. Figure[2]shows 
the LWRS spectrum over the full spectral range ~ 900-1200 A. 
Apart from the lines of neutral species, a remarkable detection 
is the interstellar O vi line in NGC604, tracing hot gas, at a ra- 
dial velocity of -211 + 25kms-'(^ = 1031 A). We did not 
observe any significant differences in the extracted spectrum 
whether correcting for the jitter of the satellite or not. 

We investigated the neutral gas inhomogeneity by using 
HST/STIS spectra toward individual stars in the ionizing clus- 
ter ("Sect. l5~4l . These spectra were taken with the G140L grat- 
ing using the FUV-MAMA detector, and have been kindly pro- 




Fig. 2. FUSE LWRS spectrum of NGC604, indicating the most prominent interstellar lines from NGC604. We also show the wavelength 
ranges of the two broad O vi and C m stellar P Cygni profiles. 
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vided to us by F. Bruhweiler and C. Miskey. The extraction 
technique is described by Miskey et al. (2003). 

We also used IUE archival data in the short-wavelength 
low -resolution mode to determine the H i content in NGC604 
using the Lya line ("Sect. 15.3V We analyzed 10 IUE spectra with 
exposure times longer than ~ 5 ksec. 



3. Data analysis 

The most plausible physical situation in NGC604 is the pres- 
ence of clouds with different physical properties that pro- 
duce different absorption line components. Indeed, the multi- 
ple sightlines toward the massive stars in the NGC604 cluster 
contribute to the global spectrum of the region. However, due 
to resolution effects, we only observe a single component, and 
this is the assumption we make in a first analysis presented in 
this section. Jenkins et al. (1986) found that this assumption, 
even in the case of complex sightlines, does not yield signif- 
icant systematic errors on the column density determination, 
unless very strongly saturated lines are present. We discuss this 
issue in detail in Sect.|9] 

3.1. The profile fitting method 

The interstellar absorption lines we observe appear symmet- 
ric and can be reproduced by a Voigt profile. We detect lines 
arising from both NGC604 and the Milky Way clouds ly- 
ing along the sightline. The NGC604 absorption component 
is blue-shifted by ~ -250 km s _1 (~ 0.8 A) with respect to the 
Galactic component (~ 0kms~'). Another weaker absorption 
component is found at ~ -150 kms -1 , which is likely to be 
associated with a high velocity cloud in M33. 

The data analysis was performed using the profile fit- 
ting procedure Owens (Lemoine et al. 2002) developed at the 
Institut d'Astrophysique de Paris by Martin Lemoine and the 
FUSE French Team. This program returns the most likely val- 
ues of many free parameters such as temperatures (T), radial 
velocities (v), turbulent velocity dispersion (b), and column 
densities (AO of species, by a x 2 minimization of the fits of 
the absorption line profiles. The program also allows changes 
in the shape and intensity of the continuum, in the line broad- 
ening, and in the zero level. The errors on the N, b, v, and T 
parameters are calculated using the Ay 1 method described in 
Hebrard et al. (2002), and they include the uncertainties on all 
the free parameters (notably the continuum level and shape). 
All the errors we report are within 2 cr. 

The Owens procedure is particularly suited to far-UV spec- 
tra overcrowded with absorption lines, since it allows for a si- 
multaneous fit of all lines in a single spectrum. This method 
allows us to investigate blended lines with minimum system- 
atic errors: for example, an Fe n line that is blended with an H2 
line can still be used to constrain the Fe 11 column density, since 
the H2 column density is well-constrained by all the other H2 
lines fitted simultaneously in the spectrum. 

We derived column densities from profile fitting using two 
different approaches: 



- Simultaneous fit. Species are arranged in several groups, 
each group defined by common b, T, and v parameters. 
One group refers to species that we assume to be mainly 
present in the neutral phase (i.e., H 1, N 1, 1, Si n, P n, Ar 1, 
and Fe n), another group corresponds to the molecular hy- 
drogen H2, another one to the interstellar O vi, and the last 
group refers to the other minor species, mainly higher ion- 
ization states. 

- Independent fits. Species may not entirely coexist in the 
same gaseous phase, hence we no longer assume similar 
b, T, and v. Each species is defined by its own physical 
parameters. 

We are able to use the NGC604 spectra to compare the 
two approaches and to discuss in particular the reliability of 
the simultaneous-fit method. This one was indeed preferred for 
spectra of blue compact dwarf galaxies (see the references in 
the introduction) - generally because of a relatively low signal- 
to-noise ratio - but could introduce systematic errors with re- 
spect to the more realistic independent-fit method. 

3.2. Stellar contamination and continuum 

In our observations, the continuum arises from the combination 
of the spectra of numerous UV bright stars. 

Contamination by stellar absorption lines could be an issue 
for the interstellar line fitting, especially for H 1. However, the 
contribution of the photospheric H 1 lines becomes significant 
only when early B stars dominate (Gonzalez et al. 1997), i.e., 
for bursts older than ~ lOMyr (Robert et al. 2003). The age of 
the burst in NGC604 is ~3-5 Myr (see introduction). Hence, 
given this relatively young age, we expect the photospheric 
H 1 lines to be relatively narrow, and their contribution to be 
negligible as compared to the interstellar damped profile. We 
verified this assumption by comparing the observed FUSE and 
HST/STIS spectra with an ad hoc Starburst 99 synthetic model 
of a young stellar population. The latter was kindly provided 
by C. Leitherer and F. Bresolin (private communication), and 
it makes use of a new library of theoretical stellar atmospheres 
that does not include interstellar absorption lines. The adopted 
spectra assume an instantaneous 3.5 Myr old burst and a gas 
metallicity of 0.4 Z . The interstellar H 1 lines (see for instance 
Ly/3 at 1025.7 A and Lya at 1215.7 A) are not significantly con- 
taminated by the - relatively narrow - photospheric lines (see 
Fig. [5}. Notice that the region around the O vi (~ 1040 A) stel- 
lar line is not reproduced well because the input physics is still 
missing from the theoretical stellar atmosphere models we use. 

This synthetic spectral model was also used as a safety 
check to compare with - but not to constrain - the continuum 
we adopted for the line profile fitting with Owens (see next sec- 
tions), and no significant differences were found. 

3.3. Influence of the apparent extent of the UV cluster 
on the line profiles 

Interstellar absorption lines are broader in the FUSE LWRS 
spectrum than in the MDRS spectrum. The line broadening re- 
sults from the convolution of several effects: 
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Fig. 3. Theoretical stellar spectra (thick solid lines) are overplotted onto the observed FUSE and HST/S77S data (histograms). The flux (in 
arbitrary units) is plotted versus the wavelength (A) for the three spectral regions between 900 A and 1200 A (FUSE wavelength range), and 
for the spectral window 1200-1800 A (STIS wavelength range). 



The intrinsic line width, due to the thermal and the turbulent 
velocities. Considering the large turbulent velocities usu- 
ally measured in BCDs and in NGC604 (see Sect.l6.1>. the 
thermal component should be negligible: 

Av = ^ x J^-+b 2 *^xb, 

" c y m c ' 

where vo is the frequency and m the mass of the species 
under consideration. 

The radial velocity dispersion of the clouds within the aper- 
ture. This effect is expected to be weaker for a single Hn 
region than for observations of entire galaxies containing 
several H n regions. 

The instrumental line spread function (LSF). For bright 
point-like sources observed with FUSE, the FWHM of the 
LSF is ~ 12 pixels, i.e., ~ 20kms _I (Hebrard et al. 2002). 
The misalignments of the individual sub-exposures. The fi- 
nal spectrum results from the co-addition of several expo- 
sures with possible wavelength shifts between individual 
spectra. This unavoidably introduces some misalignments 
because of the relatively low signa-to-noise ratio of each 
exposure. However, given the relatively good quality of the 
NGC604 observations, the misalignments do not yield a 
significant additional broadening. 



- Finally, the spatial distribution of the UV-bright stars along 
the slit dispersion axis that results in wavelength smearing. 

When comparing the FUSE spectra with the two different 
apertures, only the geometric effect due to the distribution of 
the massive stars within the apertures is likely to vary signifi- 
cantly. The Owens procedure makes it possible to estimate the 
line broadening parameter, which accounts for all the effects 
mentioned above, except the intrinsic line width, implying the 
turbulent velocity and the temperature, which is treated inde- 
pendently. Hence, we have here the opportunity of quantifying 
in a first approximation the aperture effects related to the ob- 
servation of an extended source. 

The most likely value of the line-broadening parameter for 
the MDRS (4"x20" aperture) spectrum is ~ 15 pixels (i.e., 
~ 25 km s ), which is close to the point-source LSF. By taking 
into account that a 30" source spreads the line by 100 kms 1 
(60 pixels), this translates into a spatial extent of about 5", 
which is comparable to the size of the MDRS aperture in the 
dispersion direction. The most likely broadening value for the 
LWRS (30"x30") data is instead about 30 pixels. This is equiv- 
alent to a spatial extent of 16", which is smaller than the size of 
the aperture and therefore would correspond to the actual ex- 
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Wavelength (Angstroms) 

Fig. 4. Detection of molecular hydrogen lines (at ~ 1048 .4 A 
and ~ 1052.3 A) arising in NGC604. The thick line represents 
the fit to the data (histogram, FUSE MDRS spectrum). Lines 
labelled between parentheses refer to the Milky Way compo- 
nent, the others arise in NGC604. 
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Fig. 5. Excitation diagram of the diffuse H 2 in NGC604. Column den- 
sities of the H 2 rotational levels divided by their statistical weights are 
plotted against their excitation energies. The ratio of J = 1 to J = 
levels yields a rotational temperature r=112±10K. 



tent of the UV bright stellar cluster (see the HST/S77S image at 
2000 A for comparison in Fig.^. 

As a conclusion, the line broadening parameter we infer 
for the FUSE observations is consistent with the extent of the 
source at the observed wavelength as constrained by the aper- 
ture size. The column densities derived in Sect. 16.21 take this 
broadening effect into account. 

4. Molecular hydrogen 

FUSE gives access to many molecular hydrogen lines of the 
Lyman and Werner bands. The H2 lines arising from Galactic 
clouds along the sightline are clearly present in our spectra at 
the radial velocity of 4 + 5kms _1 . These lines are responsible 
for blending other atomic lines, including lines of the neutral 
species in NGC604. However, the column density of each ro- 
tational state of the Galactic H2 is well-constrained given the 
large number of lines available, allowing correction for this 
blending effect. 

In contrast to Bluhm et al. (2003), we detected H2 in 
NGC604 in both LWRS and MDRS spectra (see Fig.0). This 
positive detection has been enhanced by the opportunity given 
by the Owens procedure to fit simultaneously - and thus gath- 
ering all the information of - all the H2 absorption lines. 

The velocity we infer for H 2 lines in NGC604 is -250.4 ± 
2.1 kms -1 in the LWRS observation and -252.2 + 3.4kms _1 
in the MDRS. We did not detect H2 from the weak absorption 
component at ~ -150kms _1 . In Table [2 we report the detec- 
tion levels and the column densities of each rotational state. We 
used these column densities to build the excitation diagram of 
Fig- HI The ratio between the H2 column densities in the 7=1 
and 7 = levels (ortho- to para-hydrogen) yields a rotational 
temperature of T — 112 + 10 K. This temperature can be iden- 
tified as the gas kinetic temperature, since these two rotational 



states should be populated mainly by collisional processes. The 
higher excitation of levels 7 = 3 and 7 = 4 is common in the 
ISM and is due to non-collisional processes, such as UV photon 
pumping, shocks, and formation of H2 on dust grains. 

Given the fact that upper rotational states (7 > 4) are not 
detected, we choose to neglect their contribution to calculate 
the total molecular hydrogen column density which is found to 
be logMFL;) = 16.86^034. The molecular fraction defined as 
/ H2 = 2 x N(U 2 )/[2 x N(H 2 ) + N(Hi)] is then 2.6 x 10~ 4 (see 
Sect.|5]for the adopted value of the H 1 column density). Such a 
low value can seem somewhat surprising in a star-forming re- 
gion considering that H2 is an important gas reservoir in which 
to form stars. Hoopes et al. (2004) found similarly low frac- 
tions in their sample of starburst galaxies, where upper limits 
of f m range from 4.9 x 10~ 6 to 1.6 x 10~ 4 . It is likely that the 
diffuse molecular hydrogen is destroyed by the incident UV 
flux from the massive stars in the star-forming region. Most of 
the remaining molecules should be in dense clumps, which are 
opaque to far-UV radiation, and do not contribute to the ob- 
served spectra (Vidal-Madjar et al. 2000; Hoopes et al. 2004). 



5. Neutral hydrogen toward NGC604 

The broad interstellar H 1 absorption measured toward NGC604 
cluster results from the blended lines arising in the Milky Way, 
in NGC604, and in the weak component at ~ -150kms _1 . 
Note that, given the intermediate velocity of the latter, its H 1 
absorption falls in the heavily saturated core of the main H 1 ab- 
sorption line. Therefore, this component does not significantly 
modify the integrated H 1 absorption line profile along the sight- 
line. 
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H2,J=0 


H2.J=1 


H2,I=2 


H2,J=3 


H2J=4 


H 2,tot 


LWRS 

detection level b 


16 39 +0 20 
~ 12a- " 


16 51 +a23 

lu.ji _ 031 

~ 14a- 


15 12 +a49 
~ 7a- 


16 17+ 030 

iu.i '_ 0S7 

~ 13a- 


14 40+ - 48 

-0.60 

~ 4a- 


16 - 86 -034 


MDRS 


16.08*$ 


16.60^ 


15.08!^ 


16 29 +a22 


14 3Q+U2 
it.jo_ 330 


16.86!°;i 


detection level b 


~ 6a- 


~ 6a- 


~ 2.5a- 


~ 6a- 


~ 2a- 





a Value obtained by summing over all detected rotational states. 
b Calculated with the Aj^ 2 method, using all the observed lines. 



5.1. Galactic component 



Interstellar clouds of the Milky Way contribute to the ab- 
sorption line spectra toward NGC604. However, according to 
Velden (1970), the Galactic Hi column density is expected to 
be somewhat lower than the intrinsic H i component. From the 
survey of Heiles (1975), the Galactic H i column density should 
be a few 10 20 cirT 2 . This is consistent with the estimates we ob- 
tained from the reddening of E(B-V)=0.09 (Israel & Kennicutt 
1980). Indeed, by assuming that jV ra /E(B-V) = 5.8 x 10 21 
cm" 2 mag 1 (Bohlin et al. 1978), we find a Galactic Hi col- 
umn density of log N(H i) m 20.7 (where the column density 
is expressed in atom cm -2 ). This value agrees well with our 
estimate obtained from the FUSE data (see Sect. 15. 2> . 



5.2. FUSE observations of Lyfi 

The Hi Lyman lines from LyyS to Ly/v fall within the FUSE 
wavelength range. However, in order to infer the H i column 
density from FUSE spectra, we can only reliably use Ly/3, since 
it is the only H i line showing damping wings. The other Lyman 
lines fall in the flat part of the curve of growth, where there is 
degeneracy between the turbulent velocity and the column den- 
sity. In addition, the higher-order H i Lyman lines are located in 
spectral regions that are more crowded by other contaminating 
absorption lines (e.g., Galactic H2 or other atomic interstellar 
lines). 

Nevertheless, the H 1 Lyman /? line in NGC604 is contam- 
inated by the stellar Ovi P Cygni doublet at 103 1.9 A and 
1037.6 A, as also observed by Cannon et al. (2004) in the 
starburst galaxy NGC625. The two stellar lines are heavily 
blended, resulting in a broad P Cygni shaped line. In an at- 
tempt to model the stellar O vi feature in order to further es- 
timate the interstellar H 1 absorption, we decided to reproduce 
the absorption component of the total P Cygni profile using a 
Gaussian profile. This is a rough approximation since the the- 
oretical asymmetrical P Cygni profile should be used instead. 
However, one can reasonably expect the global profile, which 
is the combination of the individual profiles of different types 
of stars, to have a roughly symmetrical absorption component 
(the C m, C iv, and Si iv stellar absorptions indeed appear sym- 
metrical in the FUSE and IUE spectra). We performed a si- 
multaneous fit of the interstellar lines, in particular the H 1 line 
Ly/? in NGC604, together with the broad O vi stellar absorp- 
tion. The velocity shift and the width of the latter were allowed 




1030 1040 
Wavelength (Angstroms) 



Fig. 6. The interstellar Ly/j profile as observed in the FUSE LWRS 
spectrum is heavily contaminated by the stellar O vi P Cygni doublet 
absorption. The emission component of the latter was not accounted 
for in the fitting (see text). For display purposes, all the other interstel- 
lar absorption lines are not marked, and data are rebinned by a factor 



to vary without any constraints. The emission component of the 
P Cygni profile was ignored for the fitting. 

The best result, corresponding to the minimum x 1 , is 
shown in Fig. [6] The stellar absorption has a likely width of 
~ 1300kms -1 and a blue shift of ~ -1500kms -1 . An Hi col- 
umn density of log jV(H i) = 20.7 was derived for the Milky 
Way component (estimated uncertainty 0.3 dex at 2 <x), while a 
value of logjV(Hi) = 20.75 (+0.3, see Fig.® was inferred for 
the Hi in NGC604. 



5.3. IUE observations of Lya 

The H 1 column density can also be determined from the Lya 
absorption signature in low-resolution SWP-IUE spectra (see 
Table®. We investigated Lyo- from 10 independent spectra 
(Table®. 

For the profile fitting, we used the Galactic Hi value of 
logjV(Hi) = 20.7, which is our best guess from the FUSE ob- 
servations. In each spectrum, Lya is blended with a single stel- 
lar Nv P Cygni profile (see Fig.®. To estimate this contami- 
nation, we decided to model the stellar N v absorption compo- 
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Fig. 7. Profiles for an Hi column density in NGC604 of log/V(Hi) = 20.75 (thick line), logA'(Hi) = 20.75 ± 0.30 (dashed lines), and 
logA'(Hi) = 20.75 ± 0.60 (dotted lines). Column densities of other species (including Galactic Hi) are considered as free parameters in the 
fitting routine. 



nent using the same technique as in the previous section. The 
mean width of the stellar absorption N v in the IUE spectra was 
found to be ~ 1200 + 200km s _1 and the mean velocity shift 
1600 + 150km s~ l . As a result, we inferred a mean Hi col- 
umn density in NGC604 of log N(U i) = 21.07^. 

5.4. Hi inhomogeneities revealed by STIS 

The H i column density derived from FUSE and IUE spectra 
could be different from the true H i column density, due to the 
presence of multiple unresolved absorbers along the many con- 
tributing sightlines. We investigated this issue and inquired how 
the result is affected by the inhomogeneities in the H i by an- 
alyzing individual stellar spectra from the long-slit, spatially 
resolved, HST/STIS data. A few sightlines were impossible to 
investigate because of edge effects in the MAMA detector, to- 
gether with a location of the star at the border of the slit. 

We first performed a profile fitting of Lya in the single 
stellar spectra in order to infer the column density of the H i 
for each sightline. We adjusted the profile of the stellar N v 
line using the method described in the previous sections. We 
also assumed that the Galactic H i column density is identical 
for all the sightlines, given the relatively low angular extent 
of NGC604. The results of the profile fitting are reported in 



Table 3. Hi column density determinations (in cm 2 , logarithmic 
units) derived from the IUE spectra of NGC604. 



Dataset 


Exp. Time (ksec) 


log/V(Hi) 


swp 16034 
swp 16035 


10.8 
11.3 


20.99!^ 
21 17+° :i2 


swpl9154 


18.0 


< 21.30 a 


swpl9181 


15.3 


21-02!™ 


swp24509 


15.6 


20.85!"!! 


swp04162 


4.8 


21 21 +0 - 22 


swp05682 


6.0 


21 01 +02<) 


swp06638 


4.8 


< 21.33 a 


swp07349 


6.4 


91 97+0.07 


swp24508 


5.4 


20.87!°" 


Mean column density 


21.07!"- 



a The Hi Lya absorption is severely contaminated by a terrestrial 
airglow. 



Table 0] We find variations up to 1 dex in the Hi column den- 
sity of NGC604, as compared to the uncertainties on the order 
of < 0.4 dex, suggesting inhomogeneities of the diffuse neu- 
tral gas in front of the ionizing cluster. This could be a source 
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Wavelength (Angstroms) Wavelength (Angstroms) 



Fig. 8. Example of a IUE spectrum toward NGC604. The thick line 
represents the fit to the data (histogram). Lines marked upwards are 
stellar in origin. The other lines are interstellar. The airglow emission 
around 1210 A has been removed. The emission components of the 
P Cygni profiles have not been taken into account in the profile fitting. 

Table 4. H i column density toward individual stars of NGC604 clus- 
ter, as derived from the Lya line detected with HST/S77S. 



Star" Spectral type" Flux at -1280 A b 



lO! 



117 


04 II 


4.0 


564 


09 II 


2.8 


578b 


09 la 


4.5 


675 


07H 


2.8 


690a 


05 III 


1.8 


690b 


BO lb 


3.4 


825 


05II 


3.5 


867a 


04 lab 


10.0 


867b 


04 la 


22.0 



20 
20 
20 
21 
20 
20 
20 
20 
20 



Mean column density 20 
Flux-weighted mean column density 20 



;ff(Hi) 

26 +029 

• z,u -0.40 
nc+0.16 
■ yJ -0.19 
87+0.05 

tq+0.16 
•■' -0.13 
7Q+0.36 

62 +0.26 

JO -[32 
70+ 0.25 

71+0.29 
■ ' -0.49 

•° M).26 
77+0.23 



-0.33 



a From Miskey et al. (2003). 
b In units of x!0~ 15 ergs cm -2 : 



of systematic errors when determining the total column density 
from the integrated light of the cluster. However in our case, the 
mean value of the H i column density over each single sightline 
is comparable to the mean value weighted by the star lumi- 
nosity (see Table @), which is what we actually measure in the 
global spectrum of a star cluster. 

In addition, we synthesized the global spectrum within the 
STIS 52"x2" slit, summing each one of the extracted stellar 
spectra. In this way, we simulated the spectral observation of 
a cluster as a whole, i.e., our FUSE observations. Given the 
relatively high dispersion of the G140L grating (~ 2 A), com- 
bined with the width of the 5775' slit, different positions of stars 
within the slit can lead to in a relative wavelength shift up to 
several thousands of kms -1 . We thus corrected the data for 
these shifts. This allowed us to reach an accuracy on the or- 



Fig. 9. Profile fitting of the Lyo- line from the global HST/ST/S spec- 
trum (obtained by summing the individual extracted stellar spectra). 
See Fig.|8|for a description of the plot. 



der of 1 pixel, corresponding to a velocity dispersion of less 
than ~ 150kms _1 , comparable to the wavelength smearing 
along the FUSELWRS aperture (~ 50kms _1 , Sect.EHl. From 
the Lya profile of this summed spectrum (see Fig. [5}, we ob- 
tained a column density of logyV(Hi) = 20.67*24? f° r tne 
Hi in NGC604. This estimate is consistent within the errors 
with the "simple" and "luminosity-weighted" means reported 
in Table |4] This indicates that, in our integrated spectra, we 
are not misestimating the actual column density, even though 
the absorption arises from many sightlines intersecting clouds 
with different physical properties. 



5.5. The adopted Hi column density in NGC604 

In Table |5] we report the various determinations obtained for 
the H 1 column density in NGC604. The neutral hydrogen was 
also investigated using radio observations. Israel et al. (1974) 
noticed a correlation between Hn region locations and large 
Hi complexes in M33. By measuring the Hi column density 
with the 21 cm line in emission, Dickey et al. (1993) found 
log/V(Hi) = 21.38 + 0.01 (45" x 45" beam), while Newton 
(1980) found 21.43 + 0.50 (47" x 93" beam). The direct com- 
parison between these measurements, which sample the H 1 in 
all the region, and the Hi detected in absorption in front of 
the stellar cluster is not straightforward. Hence we only no- 
tice that the H 1 quantity measured in emission is larger than 
the one measured in absorption, which is what is expected 
given the source geometry. From now on we adopt the FUSE 
value for our abundance studies in the neutral medium, i.e., 
log/V(Hi) = 20.75 + 0.30. This value, together with these er- 
rors, agrees closely with all the other determinations. 



6. Heavy elements 

We were able to identify many absorption lines in the FUSE 
spectra arising from heavy elements in the neutral ISM of 
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Table 5. Hi column density determinations using FUSE, IUE and HST/STIS data. 



Instrument Line Remarks Value 



FUSE 


Ly/5 


Corrected for the stellar vi contamination 


20.75 ± 0.30 


IUE 


Lyo- 


Mean value of the sample spectra 


21.0V™ 


HST/S77S 


Lye 


Mean value from the individual star spectra 


20.84+»i» 


HST/STIS 


Lyo- 


Flux-weighted mean value from the individual star spectra 


20.77+og 


HST/S77S 


hya 


Profile fitting of the global spectrum 


20.67^ 


Radio 


2 1 cm line 


Quantity of H i in the whole region 


~ 21.4 



Adopted column density for metal abundance derivations 20.75 ± 0.30 




Wavelength (Angstroms) 



Fig. 10. Result of the profile fitting method for the Fe n line at 
1125.4478 A. Three absorption components are detected, NGC604, 
the line of sight cloud (LOSC) at ~ -150kms~', and the Galactic 
component. 



NGC604 (see Fig. yj. From the profile fitting of these lines, 
we derived the turbulent velocities, radial velocities, and col- 
umn densities of each species. The complete list of the lines 
we analyzed is reported in Table [6] In Fig. ^3 we show the 
example of the Al 125.4 Fe n line profile fitting. 

6.1. Kinematics 

Using the results of the independent-fit method, we were able to 
compute the turbulent velocities of each species (see Table[7}. 
For comparison, the turbulence measured in the ionized gas of 
the H n region from Ho- observations is 28.3 ± 1 .2 km s~ 1 in the 
core and 20.1 ± 2.0 kms~' in the halo (Melnick 1980). Since 
turbulence does not depend on the ionic mass and, in our case, 
dominates the intrinsic line broadening (see Sect. 13. 3> . we ex- 
pected a priori to find similar values of b for all the species. 
However, we actually obtain values that are inconsistent with 
each other within the errors. This could be due to the fact that 
for most species we are only dealing with weak lines whose 
profile does not depend significantly on b. On the other hand, 
the small error bars seem to suggest a different explanation for 
this inconsistency. The absorption lines we are detecting may 



Table 7. Turbulent velocities (b) in kms~', neglecting the temper- 
ature broadening component, where errors are given at 2cr and 'IF' 
stands for independent fits (see Sect. B.H . 



Species 


LWRS, IF 


MDRS, IF 


Ni 


27.9 


+4.7 
-3.9 


20 r 48 


Oi 


26.5 


+5.3 
-5.7 


37 9+5.4 
J '- z --6.9 


Sin 


34.8 


+5.1 
-7.3 


43 - 1+ 6 6 


Pn 


19.4 


+8.3 
-3.9 


38.9 + g-« 


An 


29.7 


+ 10.5 
-11.4 


20.4^1 


Fen 


33.8 


+2.5 
-2.2 


56.5!« 


Mean value 29 ± 


13 


36 ± 14 



Table 8. Radial velocities of the neutral species absorption lines in 
km s _I , where 'IF' stands for independent fits (see Sect. l3.1> and errors 
are given at 2 cr. 



Species LWRS, IF MDRS, IF 



Ni 


-242.7^| 


-235.7™ 


Oi 


-239.8^ 


-231.0+11 


Si n 


-245.4+JI 


-245.5 


Pn 


-252.4+Ij 


-250.8**3 


An 


-242.9^ 


-246.5+3] 


Fen 


-239.8+^ 


-229.2+f 8 



Mean value -244 ± 7 -240 ± 10 



be composed of several unresolved components with various 
line parameters (i.e., width, column density, and radial veloc- 
ity). As a result, the b value derived from a single component 
analysis possibly has no physical meaning (Hobbs 1974). 

The mean b value derived from the independent fits is sim- 
ilar to the turbulent-velocity determination from the simultane- 
ous fit, i.e., 29 + 13 vs. 30.7+1* km s" 1 for the LWRS spectrum 
and 36 + 14 vs. 24.1+^kms _1 for the MDRS spectrum. This 
demonstrates that the simultaneous fit effectively averages the 
b parameter over the different atomic species. 

In addition to the turbulent velocities, we also derived 
the radial velocities from the absorption lines of the neutral 
species. Given the relatively large number of the lines we si- 
multaneously analyzed for each species, we were able to reach 
accuracies smaller than the instrumental broadening on these 
determinations. Our results (see Table[8) suggest that the neu- 
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Table 6. Analyzed lines of the metallic species in NGC604, where / is the oscillator strength, A mst the rest wavelength (from the tables of 
Morton 1991 & 2003), and /l obs the observed wavelength as obtained with the independent fit method applied to the LWRS FUSE spectrum. 



Species A mst A ob!} f comments 



Ni 


963.9903 


963.218 


0.124x10"' 


saturated, blended with Galactic H2 




964.6256 


963.853 


0.790xl0- 2 


not saturated, blended with Galactic H 2 




965.0413 


964.269 


0.386xl0" 2 


not saturated, blended with Galactic P 11 and N 1 




1134.1653 


1133.257 


0.146x10"' 


not saturated 




1134.4149 


1133.507 


0.287x10"' 


not saturated 




1134.9803 


1134.072 


0.416x10"' 


saturated, blended with Galactic N 1 


Oi 


924.9500 


924.211 


0.154xl0" 2 


not saturated, blended with Galactic H2 




929.5168 


923.778 


0.229xl0" 2 


barely saturated, blended with Galactic H 2 and Ly£ 




936.6295 


935.881 


0.365xl0" 2 


saturated 




971.7382 


970.962 


0.124x10' 


strongly saturated, blended with Lyy 




1039.2304 


1038.400 


0.907X10" 2 


strongly saturated, blended with Galactic H 2 


Sin 


1020.6989 


1019.872 


0.168x10' 


not saturated 


Pn 


963.8005 


962.997 


0.146x10' 


strongly saturated, blended with Galactic H2 




1152.8180 


1151.857 


0.236x10° 


not saturated 


An 


1048.2198 


1047.364 


0.263x10° 


not saturated 




1066.6599 


1065.789 


0.675x10' 


not saturated 


Fen 


1055.2617 


1054.457 


0.615x10-2 


not saturated 




1062.1517 


1061.342 


0.291x10-2 


barely detected 




1063.1764 


1062.366 


0.600x10' 


not saturated 




1063.9718 


1063.160 


0.475x10-2 


blended with Galactic Fe 11 




1081.8748 


1081.050 


0.126x10' 


not saturated, blended with Galactic H2 




1096.8770 


1096.040 


0.327x10"' 


not saturated 




1112.0480 


1111.200 


0.446x10-2 


not saturated 




1125.4478 


1124.589 


0.156x10-' 


not saturated 




1127.0984 


1126.218 


0.282x10-2 


not saturated 




1133.6654 


1132.801 


0.472x10-2 


not saturated 




1142.3656 


1141.494 


0.401x10-2 


not saturated 




1143.2260 


1142.354 


0.192x10' 


not saturated, blended with Galactic Fe 11 




1144.9379 


1144.065 


0.106x10° 


saturated 



tral gas is somewhat less blue-shifted than the ionized gas in 
the H 11 region, for which the radial velocity as inferred from 
Ha observations is -25610ns" 1 (Tenorio-Tagle et al. 2000). A 
similar trend in the radial velocities was found in the dwarf 
irregular galaxy NGC1705 and was attributed to a different ori- 
gin of the two gaseous phases within an outflowing superbubble 
(Heckman et al. 2001). However, there is no evidence of strong 
gas ouflows/infall in NGC604. 

6.2. Column densities 

One usually models photodissociation regions (PDRs) as a 
single-interface Hn region/molecular cloud, and the mate- 
rial in the PDRs is generally assumed to be relatively uni- 
form. However, on larger scales, the photodissociation front is 
likely to be rather clumpy, composed by small dense molec- 
ular clumps embedded in a more diffuse and mainly neutral 
medium. The PDRs are located on the surfaces of these clumps 
(Meixner & Tielens 1993). The clump distribution controls the 
far-UV photon penetration in the sense that the far-UV radia- 
tion passes between the dense clumps and is able to partially 
ionize the diffuse gas (see Fig. II 11 1. According to the above- 
mentioned sketch, the ISM in NGC604 is probably quite com- 
plex and the absorption is likely to arise in different gaseous 



: 




Diffuse neutral gas 
n-few 10cm- 3 , T-100K 




IGM 
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* 
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PDR 



Fig. 11. An illustration of the interface between the ionized gas of the 
H 11 region and the surrounding neutral medium. 

phases, i.e., diffuse ionized gas, diffuse neutral gas, or small 
dense H2 clouds. However, the densest clouds are opaque to 
the far-UV radiation and do not contribute to the total UV 
spectrum. Moreover, we cannot see the direct contribution of 
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the PDRs to the absorption spectrum, since the UV radiation 
generating these PDRs is blocked by the foreground associated 
dense molecular clouds. On the other hand, this contribution 
could manifest itself in scattered light, since dust is certainly 
present on the surfaces of the dense clumps (Mafz-Apellaniz et 
al. 2004). However, the presence of scattered light has not yet 
been assessed. We thus expect the neutral absorption lines to 
arise mainly from the diffuse neutral gas in NGC604, with a 
possible contribution from the PDRs. 

In Table |9] we report the column densities inferred from 
the two FUSE observations with both methods, as explained in 
Sect. 13. II i.e., independent fits (IF) for each element and simul- 
taneous fit (SF) for several grouped elements. Determinations 
of Sin, Pn, and Fen column densities, using the two ap- 
proaches with the same LWRS spectrum, are consistent within 
the errors, while those of N i, O r, and Ar i are significantly dif- 
ferent. The disagreement could be due to systematic errors in- 
troduced by the SF approach. This method assumes that species 
share temperatures, heliocentric velocities, and turbulent veloc- 
ities. The errors associated to this assumption are not included 
in the uncertainties we mention. 

The IF method was not used to determine column densities 
from the MDRS observation because of a signa-to-noise ratio 
too low that results in unstable solutions. Furthermore, the re- 
sults from the two different observations (LWRS and MDRS 
apertures, SF method) do not in general agree with each other 
well. The inconsistencies could be explained by the different 
aperture sizes, together with the large extent of NGC604 clus- 
ter in the far-UV (see Sect.|2}, or by the systematic effects dis- 
cussed in the previous paragraph. 

In their study of the diffuse molecular hydrogen content in 
M33, Bluhm et al. (2003) also derived Oi, An, and Fen col- 
umn densities in NGC604 (and associated lcr errors) by using 
the curve-of-growth method based on the equivalent widths in- 
ferred from the same FUSE LWRS spectrum we analyzed. The 
authors obtained 16.201"q2q for O i, which is consistent with our 
LWRS, SF value, but significantly lower than the LWRS, IF 
value. Their derived An column density, 13.65^q[q, is also 
lower than our determinations. Finally, the estimate of the Fe n 
column density, 15.00^J ][j, is only marginally consistent with 
our values. 

7. Modelling the ionization structure with CLOUDY 

In order to derive abundance ratios from column density mea- 
surements, it is generally assumed that the primary ionization 
state of one element is representative of its total abundance 
in the neutral gas. We expect to find all elements with larger 
ionization potential than that of hydrogen (13.6 eV) as neu- 
tral atoms in the Hi gas. This is the case for N, O, and Ar, al- 
though some fraction of argon (and to a lesser extent nitrogen) 
can be singly ionized in low-density neutral regions, due to a 
large photoionization cross-section (Sofia & Jenkins 1998). Fe, 
P, and Si are mostly found as single-charged ions, with negligi- 
ble amounts of neutral atoms. Thus Ni, Oi, Sin, Pn, An, and 
Fe n should be the dominant forms of the respective elements 
in the neutral gas, and their column densities are thought to be 
representative of the abundances of the element. 



However, ionization corrections may be needed since these 
atoms or ions can also exist in the ionized gas of the H n region 
along the sightlines, contributing to the absorption lines we ob- 
serve. In order to estimate this contamination, we modelled the 
ionized gas using the photoionization code CLOUDY (Ferland 
1996; Ferland et al. 1998). We assumed that the H n region is a 
homogeneous Stromgren sphere (with a radius R s ), ionized by a 
single star having the same radiation field as the stellar cluster. 
Although this is a very idealized situation, this is certainly suffi- 
cient for our purpose of obtaining rough estimates of the ioniza- 
tion corrections. The input N, O, Ar, and Fe abundances in the 
ionized gas are the observed abundances. They are taken from 
Esteban et al. (2002), which accounts for electronic tempera- 
ture fluctuations and provides consistent abundances of these 
elements in the ionized gas from the same dataset. The input P 
and Si abundances are calculated assuming, respectively, that 
P/O is equal to the solar ratio (see Lebouteiller et al. 2005 for 
P/O measurements in the Milky Way) and that Si/O is equal to 
the mean value measured in the ionized gas of BCDs (Izotov et 
al. 1999). The hydrogen volumic density is from Melnick et al. 
(1980). We used two different stellar continua to constrain our 
models. In model (1), we use a stellar continuum built upon the 
observed flux at all wavelengths. Model (2) simply assumed a 
Kurucz stellar continuum at a temperature of 48 000 K. 

The resulting optical emission-line intensities obtained 
with both models are given in Table ^3 For comparison, we 
report the values of three other giant H n regions. It can be seen 
that some lines differ significantly for each object, allowing us 
to constrain the model. Results of the two models agree rea- 
sonably well with each other and with the observed intensities. 
From these models we calculated the relative ionization frac- 
tions of each species within the Stromgren sphere (see Fig.ll2i 
and derived the expected column densities of N i, O i, Si n, P n, 
An, and Fe n in the ionized gas (assuming that half of the ma- 
terial is in front of the stellar cluster). These quantities are sub- 
tracted from the observed FUSE column densities in order to 
obtain the final column densities in the neutral gas alone (see 
column 3 of Table fTTTi . The corrections are relatively small ex- 
cept for Si ii and P n. This compares well with the findings of 
Aloisi et al. (2003) in IZw 18. The authors find that ionization 
effects are negligible for H, N, O, Ar, and Fe, the only excep- 
tion being silicon. 

Note that we could not estimate the amount of species in 
higher ionization stages in the neutral gas. Indeed An, and to 
a less extent N i, could be partly ionized into Ar n and N n in 
a low-density ISM, while hydrogen is still into Hi. For this 
reason, we might underestimate the argon and nitrogen abun- 
dances in the neutral gas when using Ar i and N i for Ar and 
N. Depending on the hardness of the radiation field, Ar/H can 
be larger by 0.2 up to 0.7 dex than An/Hi (Sofia & Jenkins 
1998). The situation is expected to be much less severe for N i, 
however. 

8. Chemical abundances in the neutral and ionized 
gas 

In Table II II we report the abundances toward NGC604 in the 
neutral gas and, for comparison, in the ionized gas of the H n 
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Table 9. Metal column densities in NGC604 from the LWRS and MDRS spectra (in logarithmic units of the column density in cm -2 ), where 
errors are given at 2 <r, 'SF' stands for the simultaneous fit method, 'IF' stands for the independent fits (see Sect. l3.lt . and 'LOSC stands for 
the line of sight cloud at 150 km s -1 . 



NGC604 
LWRS, IF 



NGC604 
LWRS, SF 



NGC604 
MDRS, SF 



NGC604 
"multi" 3 



LOSC 
LWRS, SF 



Galactic b 
LWRS, SF 



Ni 


15 31+ - 34 

iJ,J1 -0.17 


15.14 


Oi 


16 52 + ° 19 


16.26 


Sin 


15 54 + °" 

iJ - J -0 06 


15.52 


Pn 


13.70+°- 


13.66 


An 


13.86+°- 


13.95 


Fen 


14.89+°- 


14.88 



+0.02 
-0.02 
+0.05 
-0.04 
+0.04 
-0.04 
+0.06 
-0.07 
+0.04 
-0.04 
+0.02 
-0.02 



15 26+ 006 

1J,iu -0.05 

16 40+° 34 

lo ' w -0 .14 

15-62!°;i? 
13.68+°| 4 
14 06 + ° 07 
1 4 04+° 04 



[15.1, > 16.0] 
[16.3, > 17.0] 
[15.5, > 16.5] 
[13.7, 14.0] 
[13.7, 14.0] 
[14.8, 15.0] 



< 12.31 c 

14.28!°;^ 
12 75+ 032 

lz - ' J -0.58 

< 12.75 c 
14 51+ 015 



15.84 
17.13 
15.41 
14.32 
14.21 
14.95 



a Results from the multicomponent analysis are discussed in Sect.|9| 
b Errors on the Galactic estimates are on the order of ~ 0. 15 dex. 
c Upper limits calculated at 2 cr. 



Table 10. Optical emission-line intensities (normalized to /(Hy8)= 100.0) given by the observations and by our models (see text). 



Lines 


Hy 


Ha 


[Oiii] 


[Out] 


[Nil] 


[Nil] 


He i 


He i 


He i 


A (A) 


4340 


6563 


4959 


5007 


6548 


6584 


4471 


6678 


5876 


NGC604, model (1) 


47.1 


291.2 


90.3 


260.9 


11.4 


33.6 


4.1 


3.3 


11.7 


NGC604, model (2) 


47.0 


294.0 


81.8 


236.2 


10.0 


29.7 


4.2 


3.4 


12.1 


NGC604, Esteban et al. (2002) 


46.6 


291.0 


78.0 


250.0 


9.8 


26.3 


4.4 


3.7 


13.1 


NGC604, Kwitter et al. (1981) 


45.7 


281.8 


67.6 


208.9 


9.3 


28.2 


3.8 


2.7 


8.1 


NGC604, Vilchez et al. (1988) 


44.3 


286.0 


77.5 


207.7 


12.4 


33.5 


3.7 


2.7 


11.5 


NGC5461 


46.5 


291.0 


112.0 


352.0 


10.8 


31.2 


4.4 


3.6 


12.7 


NGC5471 


47.6 


278.0 


209.0 


640.0 


1.9 


6.2 


4.1 


2.9 


11.8 


NGC2363 


47.3 


278.0 


244.0 


729.0 


0.4 


1.5 


4.1 


2.9 


12.3 



1.5 



1 / 1 

/ Si 11 
/Fell 

/pn 
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Fig. 12. Ionization fraction of each species relative to H i, as a func- 
tion of the distance to the ionizing source. O i, N i, and to some extent 
An are well-coupled with Hi in the neutral gas (i.e., for distances 
larger than R s = 8.20 x 10 20 cm, which is the radius of the Stromgren 
sphere). 



region. For the latter we used abundances given in Esteban et 
al. (2002), instead of the older values from Kwitter et al. (1981) 



and Vilchez et al. (1988). This, because Esteban et al. homoge- 
neously derived the N, O, Ar, and Fe abundances from a higher 
resolution optical spectrum, accounting for electronic temper- 
ature fluctuations. We estimate the Si abundance in the ion- 
ized gas by assuming that the abundance ratio Si/O is equal 
to the mean value measured in BCDs (Izotov et al. 1999). 
Abundances are normalized to the solar values from Asplund 
et al. (2004) as [X/H] = log (X/H) - log (X/H) . Figure[l3] 
is a graphic representation of these data. 

In the neutral gas the abundances of all heavy elements that 
we obtain are consistent within 2<x with the oxygen abundance, 
which is on the order of 1/10 solar (see Table ll 11 1. If real (i.e., 
not driven by the large uncertainties in the H i column density), 
this consistency would be somewhat surprising, since we ex- 
pect at least Si and Fe to be depleted onto dust grains. Indeed, 
the medium we are considering should be comparable to the 
diffuse neutral medium in front of stars like £ Ophiuchi or // 
Colombae in our own Galaxy. N, O, P, and Ar are not de- 
pleted or only a little (within a factor less than ~ 3) toward 
these Galactic sightlines, while Si and Fe show depletions by a 
factor of at least ~ 3 and 10, respectively (Savage & Sembach 
1996; Snow & Witt 1996; Howk, Savage & Fabian 1999). In a 
similar way, the heavy element abundances in the ionized gas 
of the H ii region are also consistent within the errors with the 
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Table 11. Abundances of N, O, Si, P, Ar, and Fe in the neutral gas, where IC stands for ionization correction (see Sect.Q and errors are at 2a. 





Ion 


IC (dex) 


log (X/H) a 


[X/H] 


[X/H] b HII 


N 


Ni 


-0.09 


_c 7f)+0.30 

[-5.7, > -4.8] 


-1 48+ 031 
i/fO -0.31 

[-1.5, > -0.6] c 


U - JZ -0.32 


O 


Oi 


-0.08 


-4 57+ - 31 

t.j ; _ 31 

[-4.5, > -3.8] 


i 9Q+0.32 
-0.32 

[-1.2, > -0.5] c 


+0.00 + °; 2 J 


Si 


Sin 


-0.34 


-5 57 +031 

J ' J -0.31 

[-5.6, > -4.6] 


-1.08!- 
[-1.2, > -0.2] c 




P 


Pn 


-0.15 


-7 74+U-31 
' - x -0.32 

[-7.2, > -6.9] 


-0 60 +u 33 

u - u -0.32 

[-0.6, -0.3] c 


/ d 


Ar 


Ari 


-0.05 


-6.85+JJil 
[-7.1, > -6.8] 


-1 03 +0M 

1 '"-'-0.33 

[-1.3,-1.0]° 


+0 25 +0 - 20 


Fe 


Fen 


-0.04 


t q i +0.30 

-0.30 

[-6.0, > -5.8] 


-1 36 +0 - 31 

l,D -0.31 

[-1.4,-1.2]° 


-1 02 +0 - 22 

1-u -0.22 



a Values derived from the simultaneous fit method using the FUSEUWRS spectrum (see Sect. l6.2> . after respective ionization correction. 

b Ionized gas abundances from Esteban et al. (2002). 

c Values referring to the multicomponent analysis presented in Tableland discussed in Sect. [5] 

d No direct abundance determinations exist in the ionized gas. 
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Fig. 13. Abundances of N, O, Ar, and Fe in the neutral and ionized 
gas, compared to the solar values (from Asplund et al. 2004). We use 
the notation [X/H] = log (X/H) - log (X/H) . Empty circles are values 
for the ionized gas from Esteban et al. (2002). Filled circles indicate 
values in the neutral gas as determined from the simultaneous fit in 
the FUSE LWRS spectrum, using a single-component analysis. They 
may be underestimated for N, O, and Ar (see Sect. 9). 



O abundance, which is solar in this case. The only exception is 
the gaseous Fe, which is about 1/10 solar. 

If real, the observed abundance trends indicate that nitro- 
gen, oxygen, and argon are lower by > 1 dex in the neutral gas 
phase as compared to the ionized one. This confirms a simi- 
lar trend already observed in blue compact dwarf galaxies, and 
it contrasts with the most recent finding of no offset in the a- 
element content of the neutral and ionized gas in the nearby 
damped Lyman-a galaxy SBS 1543+593 (Bowen et al. 2005). 
Ionization corrections in the neutral medium cannot fully ex- 



Fig. 14. Abundance ratios in the neutral and ionized gas of NGC604. 
See Fig.^|for a description of the plot. The abundances in the neutral 
gas weres determined using a single-component model. 



plain the similar offsets of N, O, and Ar in NGC604, since Ar 
is expected to be much more affected than N and O (Sect.0. 

Iron is instead the only element showing similar abun- 
dances in the two gaseous phases. The iron behavior requires 
that this element has the same gas abundance and depletion 
factor in both the neutral and ionized gas phases. This is rather 
surprising, because we expect iron to be underabundant in the 
H i region and depletion on grains to be similar or larger in the 
neutral medium compared to the H n region (but see Sect. 9). 

For what concerns the relative abundance ratios, we ob- 
serve the following behaviors (see Fig. 1141 . The ff-elements O, 
Si, and Ar are roughly in solar proportion in both the neutral 
and the ionized ISM. This seems to indicate that Ar does not 
suffer from the ionization effects typical of a low-density neu- 
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tral medium (see Sect. and that Si is not more depleted in 
the H i gas than in the H n gas, another rather surprising result. 
In a similar way, the consistency of the (subsolar) abundance 
ratio N/O in both the neutral and ionized gas suggests that N 
is affected neither by ionization nor by depletion effects. The 
behavior of phosphorus in the neutral gas is also particularly 
intriguing, since P/O and P/Si in this gas seem to be super- 
solar, in disagreement with what found by Lebouteiller et al. 
(2005): these authors derived essentially solar P/O ratios for 
several sightlines sampling the Galactic diffuse ISM, as well as 
for a few damped Lyman-ff systems and the ISM of the small 
Magellanic cloud toward the star Ski 08. 

We show in the next section that the N i, O i, and Si n col- 
umn density determinations may be severely affected by sat- 
uration of hidden velocity components. If this is indeed the 
case, the analysis of Sect. 6, which assumes a single absorbing 
neutral component with a large velocity dispersion in front of 
NGC604, would result in an underestimate of the abundances 
of these elements. Once hidden saturation is properly taken into 
account, it is possible to easily reconcile all the puzzling results 
from the behavior of the various heavy elements in a very sim- 
ple way, i.e., by considering the possibility that the neutral and 
ionized gas could have similar abundances. 

9. Multicomponent analysis 

Instead of arising from a single neutral cloud with a large veloc- 
ity dispersion, the interstellar lines we detect might be the blend 
of many unresolved absorption components whose widths and 
velocities are related to the ISM structure and to the source 
morphology. As we now show, this may result in a severe un- 
derestimate of the column densities for the most saturated lines 
if a single-component analysis is made. It is not possible with 
the present data to assess whether multiple saturated compo- 
nents contribute to the absorption lines we observe. However, 
it is possible to perform empirical tests and to calculate the un- 
certainties associated with the presence of multiple - possibly 
saturated - components. A detailed analysis can be found in 
Lebouteiller (2005). In this section, we summarize the method 
and the results. 

9.1. Presence of hidden saturated components 

To identify which lines can suffer from hidden saturation, we 
first consider the presence of a test component placed at various 
velocities within the broad line. We choose for this component 
a b- value of 2 km s _1 , which we consider as the minimum value 
allowed for a single sightline intersecting a single interstellar 
cloud (Tumlinson et al. 2002 measured b values systematically 
larger than 1.6kms _1 toward stars of the Magellanic Clouds). 

Among the detected lines in NGC604, Fen Al 142 is the 
only one for which the test component can never be satu- 
rated due to the relatively low value of the oscillator strength. 
The other lines of Fe n and all the lines of the other species 
could suffer from hidden saturation by a component with b = 
2 km s . Hence, the column density derived from the observed 
global Al 142 Fen line should not be affected by the system- 
atic errors associated with hidden saturated components, even 



in the center of the global line, where the optical depth is max- 
imum. Note that it can still suffer from uncertainties due to the 
non-linear sum of the individual spectra (see previous section). 
The Fen column density obtained using the Al 142 line alone 
is logA^Fen) = 14.97 + 0.07, slightly higher than the value 
we found earlier using all the available lines (14.88*q'q2, see 
Sect. 16.21 . This confirms that using stronger Fe n lines in ad- 
dition to /U142 introduces systematic errors if there are hid- 
den saturated components. However, the errors do not exceed 
~ 0. 1 dex in the case of Fe n. 

For other species, we cannot rule out that even their weakest 
observable lines could be made of saturated individual compo- 
nents. Although we have no information on the actual velocity 
component distribution responsible for the global profile, it is 
still possible to test various plausible distributions, which can 
somehow be constrained if we simultaneously adjust and fit all 
the lines of a given species, from the weakest line to the heavily 
saturated one. 

9.2. Test of possible component distributions 

We consider here various component distributions, fully de- 
fined by the number of components n c , by the turbulent veloc- 
ity b of each component, and by the spacing AV between them 
(which can possibly be lower than b, depending on the number 
of components). 

The ideal species to test for plausible distributions is Fe n 
because of the large number of lines available, with a wide 
range of oscillator strengths. Fen might have a space distri- 
bution differing slightly from that of the other species, due to 
abundance and ionization inhomogeneities. However, here we 
only wish to test the method, so we ignore this problem. 

We consider various distributions with n c varying between 
1 and 20 and b between 2 and 44 km s _1 (which is the highest 
b value for a single component, see Sect. 16. 1> . while AV is set 
so that the components are uniformly distributed within the ob- 
served Fen line widths. For a given line, the column densities 
in each component are considered as free parameters by the 
fitting procedure Owens. 

We find that many component distributions provide a x 2 
equal to or lower than the x 1 f° r a single-absorption compo- 
nent. We thus consider all these distributions as mathemat- 
ically and physically plausible. The distributions give a to- 
tal column density (sum of all the components) in agreement 
within 0. 1 dex with the determination using a single compo- 
nent. Furthermore, we notice that satisfactory distributions im- 
plying components with low b values (potentially responsible 
for saturation effects) also imply a large number of components 
n c , so that the total column density is dispersed in many com- 
ponents with relatively low column densities. Typically, at least 
13 components with b = 2km s were needed to adjust the 
lines. We do not find satisfactory distributions implying high 
column densities together with low ^-values. 

We used the velocity distributions that we found satisfac- 
tory to fit the Fe n lines, to also fit the lines of the other species 
(Ni, Oi, Sin, Pn, and Ari). We discovered that the distribu- 
tions do not introduce strongly saturated components in the P n 
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and Ar i lines, even when b — 2 km s . The column density we 
determined using the velocity distributions is similar to the es- 
timate using a single-component fit within ~ 0.2 dex. However, 
several possible distributions imply that even the weakest N i, 
O i, and Si n lines are made of saturated components. For these 
species, the sum of the column densities in the individual com- 
ponents can be larger than 2-3 dex in comparison with the 
result of the single-component analysis. 

Thus the column densities of Ni, Oi, and Sin could be 
severely underestimated. On the other hand, the Fe n and, to 
a lesser extent, the Pn and An lines do not seem to be strongly 
affected by systematic errors due to the presence of hidden sat- 
urated components. We show the ranges of column densities 
derived from the multicomponent analysis in Table |9] 

Focussing on the P n, Ar i, and Fe n column densities only, 
the interpretation of the abundances in the neutral gas dif- 
fers from what is discussed in Sect. [8] Iron has similar abun- 
dances in the ionized and neutral phases, suggesting that the 
two gaseous phases could indeed have similar metallicities. 
The underabundance of Ari could then be explained by the 
fact that argon can be partly ionized into Ar n in this gaseous 
phase. The actual Ar/H would be closer to Arn/Hi. This re- 
sult would agree with the findings in IZw36 (Lebouteiller et al. 
2004). Finally, little can be said about the abundances of Ni, 
O i, and Si n, which could well be similar to those in the H n re- 
gion. To conclude, the evidence of a difference in abundances 
in the H i and in the H n gases essentially vanishes as a result of 
our tests. 

10. Conclusion 

This study provides the first detailed analysis of interstellar 
lines of Hi, Ni, Oi, Sin, Pn, An, and Fen in the neutral 
medium in front of a giant H n region in the spiral galaxy M33. 
Since NGC604 is a nearby system, we have been able to per- 
form the necessary critical tests for analyzing possible selection 
effects. 

- In the frame of a simple model, we have derived column 
densities of metals in the neutral gas of NGC604 from both 
FUSE LWRS and MDRS spectra. These independent esti- 
mates allowed us to quantify the effects of the source extent 
on the interstellar absorption line profiles. 

- The continuum used for the profile fitting was checked by 
comparison with a theoretical stellar model. Besides, no 
significant contamination from stellar photospheric lines 
was found for the H i absorption lines we investigated. 

- A particular attention was given to the neutral hydro- 
gen column-density determination. The Ly/3 line from the 
FUSE LWRS and MDRS spectra is contaminated by the 
Ovi P Cygni doublet. We modelled this contamination 
to obtain the Hi column density. Also, profile fitting of 
the Lya absorption in HST/STIS spectra toward individ- 
ual stars in the cluster NGC604 reveals inhomogeneities 
in the neutral gas. We finally adopted the FUSE value of 
log Af(Hi)=20.75 with a conservative error of +0.3 dex to 
account for all the possible uncertainties. Within the errors, 
this H i column density is consistent with the IUE determi- 



nation using Lya and with the 21 cm line radio observa- 
tions. 

- By modelling the ionization structure of the Hn gas with 
the photoionization code CLOUDY, we have shown that N i, 
O i, An, and Fe n are reliable tracers of the neutral gas, in 
contrast with Si n and P n, which require ionization correc- 
tions to obtain final abundances in the neutral phase. 

Adjusting absorption lines with a single component, we 
find that N, O, Ar, and Si are underabundant in the neutral gas 
as compared to the ionized gas by factors > 10, while Fe/H is 
similar in the two gaseous phases. This result is rather puzzling, 
since iron is expected to be equally or more depleted in grains 
in the neutral gas compared to the ionized one, and there is no 
reason it should be relatively more abundant in the neutral gas. 

This led us to investigate in detail the influence of individ- 
ual unresolved components in the analysis of absorption lines. 
Using the method and results of Lebouteiller (2005), we argue 
that N i, O i, and Si n column densities can be severely under- 
estimated if there are saturated hidden components, while Fe n 
and, to some extent, Ar i and P n should be more reliably deter- 
mined. Since Fe is the only element to show similar abundances 
in the neutral and ionized gas, it is possible that all elements in- 
deed have similar abundances in both media. The underabun- 
dance of Ar would then be due to the fact that we used Ar i/H i 
to estimate Ar/H, when ( Ar 1+ Ar n)/H i should be used instead. 
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